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I. INTRODUCTION 

The Standard Model (SM) in high energy physics [l| has been remarkably successful in: describing the properties 
of elementary particles, predicting the existence of the quarks c, t and b, and the third generation of leptons r- v Tl 
the existence of the eight gluons, and the weak bosons W^, Z° before their discovery, predicting parity violating 
neutral- weak-currents, and in being consistent with all the experimental results [HQ. However, the SM falls short of 
being a complete theory of the fundamental interactions because of its lack of explanation of the probable unification 
of the fundamental interactions, the pattern and disparity of the particle masses (mass hierarchy), the origin of the 
CP violation in nature, the matter-antimatter asymmetry, the pattern of quark mixing, lepton mixing and the reason 
why there are 3 generations. 

As a partial solution to confront these deficiencies, a large number of parameters must be put in "by hand" into 
the theory (rather than being derived from first principles), such as the three gauge couplings (g%, gi and 33), nine 
fcrmionic masses (six quarks and three leptons), the Weinberg angle (0 W ), four quark-mixing parameters (CKM) and 
two more parameters in relation to the Higgs potential (p and A). 

One of the most subtle aspects of the model is associated with the Higgs sector [H . The Higgs field and its non- 
vanishing vacuum expectation value (vev) is the essential ingredient to carry out the spontaneous symmetry breaking 
(SSB) required to transform the hypothetical massless particles in the Lagrangian into the actual massive physical 
particles. However, the Higgs particle has not yet been discovered. 

In this paper we study the extension of the SM with two Higgs doublets (2HDM) that presents the challenge that 
the quartic interactions between the scalar doublets are not theoretically determined. This model is studied mainly 
for three reasons. The first one is that the 2HDM has a much richer Higgs spectrum (3 neutral and 2 charged Higgses) 
and a different high energy behavior. This makes that a lower mass than in the SM Higgs is permitted. Another 
reason may be that a different pattern of hierarchy of the Yukawa couplings is possible, because of the presence of 
two independent vacuum expectation values of the Higgs fields 5]. The third reason is that the Higgs sector of the 
Minimal Supersymmetric Standard Model (MSSM) contains two Higgs doublets, so the Higgs sectors of the MSSM 
and the 2HDM are similar and the study of the 2HDM model may give important information on the properties of 
the Higgs sector in the MSSM 0. 

In Section [II] we introduce the potential for the 2HDM in a special parametrization, and briefly discuss the SSB. 
In Sections IIIII and IIVI we present the Higgs mass matrix and its diagonalization method, the mass spectrum, mass 
eigestates, and special cases of mixing, where the Higgs masses are simply related to the parameters of the potential. 
In Section [V] we obtain and classify the constrains for the quartic couplings derived from the mass formulas, from the 
vacuum stability principle through the Lagrange multipliers method, and by imposing extreme stability conditions. In 
Section VI we numerically solve the set of the renormalization group equations from which through triviality principle 
the physical bounds of the model are determined under different conditions. Finally, Section IVTI1 is devoted to the 
presentation of the results and the conclusions. 



II. THE TWO-HIGGS DOUBLET MODEL 

In the SM the fermion masses arise, after the SSB, from the couplings between the fermions and a single Higgs 
doublet. The mass ratio of the b and t quark is of the order of 1/40. To understand in a natural way the origin of 
this difference in the values of the masses of the third generation of quarks, one can assume the existence of a second 
Higgs-doublet in the Higgs sector of the SM. In this context one assumes that the quark t obtains its mass through 
the $1 doublet and the quark b from another doublet $2 Q- In this way one can explain in a more natural way the 
hierarchy problem of the Yukawa couplings, as long as the free parameters of the new model acquire the appropriate 
values. 

The Higgs sector of the 2HDM consists of two identical (hypercharge-one) scalar doublets $1 and <I>2. There are 
several proposals for the Higgs potential to describe the physical reality in the framework of the 2HDM [|| . The 
potential we consider in this paper is compatible with Ref. [10]. It is such that the CP symmetry (charge-conjugation 
and parity) is conserved, the neutral-Higgs-mediated flavor-changing neutral currents (FCNC) are suppressed in 
the leptonic sector, and in the quark-sector they are also forbidden by the GIM mechanism [ljj in the one loop 
approximation. In the Lagrangian C in which we leave out the leptonic terms, 

C = C gf + C Km + C Y -V (1) 

the Cgf and Cain correspond to kinetic parts of quarks and bosons and they contain the covariant derivatives that 
provide the interactions among the gauge bosons and the Higgs bosons. They also give rise, after the SSB, to the 
masses of the gauge bosons (mediators of the electroweak interactions) . The fermion masses are generated from the 
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Yukawa couplings in Cy 

£ y = J2 (?ffi Li *i«Ri + g^Lfrdv) , (2) 

i-j 

between the Higgs bosons and the quarks. In Cy, the g^' d ^ are the Yukawa coupling matrices. The superscripts 
(u, d) refer to the up and down sectors of quarks, respectively and the subscripts (L, R) correspond to the left handed 
doublets and right handed singlets in the quark sector. In this paper, we will focus our attention on the potential V. 



The Higgs potential 



The Higgs potential depends on seven real parameters A*2 an d A, (a = 1..., 5) from which the five Higgs masses 
come up after the SSB. The most general renormalizable SU(2) x U(l) invariant Higgs potential, that preserves a CP 
and a Z 2 symmetry ($i — > $i, $ 2 — > —$2) is given by 



V = + A*2^*2 + Ai ($1 $1) 2 + A 2 (*J$ 2 ) 2 + A 3 ($1 $1) ($3, $2) 

+ A 4 ($1*2) ($ 2 $i) + ^A 5 ($1$ 2 ) 2 + 

For the sake of simplicity a special basis is introduced 
In this basis 

V" = + /i^B + XiA 2 + A 2 B 2 + A 3 AB + UC'D' + ^A 5 [C 2 + D' 2 ] . 
The two Higgs doublets can be represented by eight real fields fa, i = 1, ... ,8, 

1 ! fa + ifa 1 ' - 



55 + ifa 
67 + ifa 



(3) 



(4) 



(5) 



(6) 



If charge is conserved and there is no CP violation in the Higgs sector, after the SSB, the non-vanishing vacuum 
expectation values of the fields fa and fa are real, 



(7) 



In terms of the fields 



(fa) = (fa) = (fa) = 0, 
the hermitian basis is given by 



5) = (fa) 



3> =0. 



A = <t>i + 44 + <t4 + <ft, B = <pi- 

C = fafa + ifafa - ifafa + fafa 
D' = fafa + ifafa - ifafa + fafa 



■ fafa + ifafa ~ ifafa + fafa, 

■ fafa + ifafa ~ ifafa + fafa- 



and after the SSB they become 



(8) 



(9) 



(A) 



(B) 



^2; 



(C) = (D 1 ) = -v lV2 . 



(10) 



III. THE MASS MATRIX 



The conditions for the minimum of the potential are obtained from the vanishing of the first derivatives at the 



minimum 



dV 



= 0, with the condition that the matrix of the second derivatives at the minimum: 



is 



4 



positive definite. Therefore 
dV 



» 2 ®A + ^l^-B + 2X x A-^—A + 2X 2 B^-B + X 3 B^- + X 3 A^ 



<O|0,|O> 



X.D'^ + X.C'^ + X, 

dpi o<pi 



d d 



from which after some simplifications two non trivial equations are obtained 

[if + X x vf + 2X T v\ = or v x =0, (4 + X 2 v\ + 1X T v\ = or v 2 = 0, 

where 

At = (A3 + A 4 + A 5 ) . 
The mass matrix elements are obtained from the equation 

1 d 2 V 



|0) = (11) 



(12) 



M 2 - 

13 2 dfadfa 



(13) 



(14) 



and the explicit form of the matrix of the second derivatives reads 

d 2 V , 2 x d 2 A , „ , , „. d 2 B OAdA , dB dB 

^ f OA dB dB dA\ (dC'dC dD' dD'\ ( dC 3D' 3D' dC 

3 \ d4>j dfa 3<pj dfaj 5 V d4> 3 dfa d(f>j dfa J 4 \ dcfij dfa Qfa dfa 

ri 2 c" 8 2 ry 

+ (X 4 D' + X 5 C) + (A 4 C + X 5 D') 



(15) 



Using Eqs. (|12[) the 16 non vanishing matrix elements are 



1 



Mf 1 = M 2 2 = --(X i + X 5 )vi M£ 3 = 2\ 1 v(, M, 
1 



M\ x 


= M 2 2 


m! 5 


= M 2 6 




— M 2 X 


Mir 


= M? 3 



f 2 

44 



-As?;. 



-| (A 4 + A 5 ) vl M 2 7 = 2X 2 v 2 , M 2 S = -X 5 v 2 , 



M( b = M£ x = Mi 6 = M£ 2 = - (A 4 + A 5 ) v lV2 , 

(A 3 + A 4 + A 5 ) viv 2 , M 2 S = M 8 2 4 = X 5 v x v 2 



(16) 



Diagonalization of the mass matrix 



The Higgs masses and the Higgs mass-eigenstates are obtained after a suitable diagonalization of the matrix in 
Eq. ([14")) . The diagonalization of the matrix whose elements are given in Eq. (fTl)|) is performed in two steps. 

A block ordering is performed and a diagonalization of each block is carried out. The blocks are obtained by means 
of the application of two consecutive unitary transformations U\ — U\ and U 2 — U 2 



(M 2 ) B ^U 2 U 1 Mf J Ulul 

The non vanishing matrix elements of the unitary transformations are 

(E7i) n = (U x ) 25 = (U x ) 33 = (U X ) M = (U x ) 52 = (U x ) 66 = (U x ) 77 = (U X ) 8S = 1, 
(U 2 )n = (U 2 ) 22 = (U 2 ) 33 = (U 2 ) 47 = (U 2 ) 55 = (U 2 ) 6e = (U 2 ) 74 = (U 2 ) S8 = 1. 



(17) 
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After carrying out both transformations, the 8x8 matrix becomes 



' ± v± 1 i 




o 


o 


o 


o 


o 


o 


\ 


Mil 


Ml 

00 


























Ml 3 


Mf 7 






















M| 3 


M| 7 




























Mj 2 


Mi 6 






















Mh 


Ml 6 




























Mh 


Ml 




Vo 

















ML 


Mis 





(18) 



The matrix in Eq. (|18[) . is ready to easily perform the total diagonalization 

\ " ^ Ml M| 5 J ® Af 7 2 3 M 7 2 7 J ® 1, M| 2 M| 6 J ® M| 4 M| 8 J 

The next step is to perform the diagonalization of each of the submatrices in Eq. (I19[) . 



IV. HIGGS MASS-EIGENSTATES BASIS 



Let us now proceed to relate the gauge states with the mass eigenstates. 
The scalar fields in Eq. © can be represented as 



where 



<j>t = 01 + i<t>2, 4>t = 05 +#6, 01 = 03 + *04, 02 = 07 + «08, (21) 

and 

03 = — ^+/ll, 04=?/l, 07 = "^+/l2, 08=??2- (22) 

Now, the physical fields (mass eigenstates) and the Goldstones (massless eigenstates) are obtained from the gauge 

eigenstates by a unitary transformation that diagonalizes the corresponding submatrices Eq. (|T9]) . in the following 
way 

(£WJn- »■ (23 » 



where 

cos a sin a 



/ cosa sma \ Tr + T r T /„,, 



and t/^ = C/ 7 have the same form as U a . a is the mixing angle between the neutral states and 2 j i- e -i 03 an d 07 
, the /3 angle is the one between the charged states, and 7 is related with the CP-Odd states </>4 and 0g. 



, (A3 + A 4 + A 5 )fii;2 7T 7T 

tana = , y — t&n2a— — 7- — 5 ; — ,7 — , — — < a < — , (25) 

1 + (A!« 2 -A 2 « 2 2 ) 2 2' 



and 



tan/3 v 2 = (v\ + v\) , < /3 < J, 7 = /?. (26) 

Wi 2 

The resulting physical particles in the Higgs-sector are: two CP-even-neutral Higgs scalars (H°,h°), one CP-odd 
neutral Higgs scalar (A ), two charged Higgs bosons (H^), and three Goldstone-bosons (G^G ) that contribute to 
the mass generation of the gauge vector bosons W ± and Z°, respectively. 
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The mass formulas 



After the complete diagonalization, we obtain the following relations: 
1. The mass eigenvalues for (H°,h°) are 



M H°,h° = A i"i + A 2«2 ± V - X 2 v 2 f + {viv 2 X T f > 0, (27) 

2. The eigenvalues for the mass eigenstates H ± and G ± are 

Af£ ± =0, M 2 h± = -I(A 4 + A 5 )v 2 >0, (28) 

3. Finally, the mass eigenvalues for G° and A are 

M 2 = 0, M| = -X 5 v 2 > 0, (29) 

As expected, after the electroweak symmetry breaking (EWSB) the eight components of the two complex isodoublet 
fields are transformed into: two charged Higgs bosons H ± , three neutral Higgs bosons H°, h°, A a , and three massless 
Goldstone fields G°, G ± (which are transformed into the longitudinal components of the gauge bosons and Z°). 
At this level, the values of M^o and M H ± are not related to the parameters Ai, A2 and A3. This means that, 
apparently, there is a complete independence between the A , and the h°, H° , which is not all true. 
As in the Standard Model, the values of the quartic couplings are not fixed by the model. To proceed as in the 
SM [l2T |. to determine the Higgs masses, one has to consider two important physical principles. The vacuum stability 
constrains the values for the quartic couplings. To have a complete view, we invert the former equations to express 
the quartic parameters in terms of the masses of the Higgs fields. 

Ai = — -2 (Mjj Q cos 2 a + Mia sin 2 a) , A 2 = i { m h° sin2 a + M h° cos2 a ) ( 30 ) 

^1 2 



(M 2 H0 -M 2 h0 ) ^ ^ t o M 2 H± ^ _ M\ - 2M% ± x M 2 
2viv 2 



X 3 = V " H " ^ sin2a + 2^, A 4 = A 5 = -^f , (31) 



At = (A 3 + A 4 + A 5 ) = ( M *° M ' o) sin 2a. (32) 

2viv 2 



Particular cases 



To obtain Eq. (|3T))) - ([3"2l we have considered the following relations: Since Eq. (f2~T|) is equivalent to 



X lV \ + X 2 v\ ± (Aiw 2 - X 2 v 2 2 ) Vl + 2/ 2 l (33) 



and 



we obtain from Eq. (1231) 



1 



as well as 



and 



Vl + V 2 = VI + (tan 2a) 2 = — , (34) 

v cos 2a 



. (A 3 + A 4 + A 5 ) v x v 2 

yj (\iv 2 - X 2 v 2 f + (A 3 + A 4 + A 5 ) 2 ( Vl v 2 ) 2 

Xxv 2 - X 2 v 2 

cos 2a = (36) 
sj (X lV 2 - X 2 v 2 ) 2 + (A 3 + A 4 + A 5 ) 2 (v lV2 ) 2 

(M| -yc O s2a = (A ^_ A ^ 2)i ^±^ =M + A2 , f] . (37) 
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With these equations it is easy to obtain Eq. ([30 



a.- In the case when the mixing angle is a = 0, i.e., At = 0, A3 > 0, 



^=2^°' A3=2 (^ 



Ai = tt M ho. A 2 = -^A^ , A 3 =2 —2- , (38) 



M A a\ 2 n fM H ±\ 2 ( Mao X 2 



A < = b^J - 2 l^J ' Xs = "l^rJ ■ (39) 

and 

M^o - Mlo = 2 (Aiw 2 - A2U2) ■ (40) 

The parameters fix, [12 in the potential Eq. are related to the neutral Higgs particles in a very simple way, similar 
to the one between the parameters A and \x in the SM. 

- 2n\ = Mfjo = 2Aiw 2 , -2^ 2 = Af 2 = 2A 2 ^. (41) 

and the vevs satisfy 

v 2 = (vl + vl) = - + £) = _L_ ( Al M* + X 2 M% ) . (42) 
In this particular case, each Higgs particle is associated with a specific parameter Ai, A2, A3, A5. 

M H o = wiV2AT, M h o = W2V2A2", M H ± = -^=Vh, M A o = vy/m. (43) 

V2 

The degeneracy in the masses M H o, M h o implies that Xivf — A 2 tj| = 0. 
b.- In the case when the mixing angle is a — tt/2, At = 



A i= — ,U,~, X, = —j\lf r . A.; = 2 I —ILL I , (44) 



2 -^M h 2 o, A2-^M|o, A 3 -2^ ± 



M A o\ 2 n fM H ±\ 2 x (M A ^ 2 



A ' = (^J ' As= v «7 ■ <45) 

The parameters ^1,^2 in Eq. (|12[) become: 

- 2/^ 2 = M 2 = 2Aiw 2 , -2^ = M% = 2X 2 vj. (46) 

and 

v 2 = (AiA^o + A 2 Af 2 ) . (47) 

As in the former case, each Higgs particle is associated with a specific parameter Ai, and (-ff , h°) interchange places. 

M h o = v X y/ZX u M H o = v 2 y/2X 2 , M h ± = ^= M A o = uyW (48) 

In this section, we have considered the main features of various special cases for the parameter a where a decoupling 
of the Higgs bosons take place, and the case where the masses of the CP-even neutral particles coincide. 

V. VACUUM STABILITY CONSTRAINS 

Bounds due to the positive mass-values 

Due to the fact that the masses are positive, from the previous results, one gets information for the allowed values 
of the Ai parameters. 

Ai>0, A 2 >0, (A 4 + A 5 )<0, A 5 <0, A 4 <|A 5 |. (49) 
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and Eq. ([27]) implies that 



AiA 2 > I (A 3 + A 4 + A 5 ) 2 . (50) 



In terms of the masses, the conditions in Eq. (|49p and Eq. ([50)1 become trivial 

Al = 2^f ^ M "° COs2 a + Shl2 ^ > °' ^ 

A 2 = 7T2 [ M h« sin2 a + M h« cos2 «] > °. ( 52 ) 
2v 2 

AiA 2 - i (A 3 + A 4 + A 5 ) 2 = -^M% Ml > 0. (53) 

To improve previous information about the allowed values of the quartic couplings and therefore for the masses, we 
have explored the consecuences of considering the vacuum stability conditions (VSC), through the method of the 
Lagrange multipliers. 

Lagrangian multipliers method and the VSC 

Considering one restriction: Let us introduce the variables Xi and the parameters hi, defined as [l3j 

£l = |$l| 2 , X 2 = |$ 2 | 2 , X 3 = i (*J$a + *£$l) , ^ = 1 ($1*2 -$ 2 $!) , (54) 

6n = Ai, & 22 = A 2 , 6 3 3 = (A 4 + A 5 ), 644 = (A 4 -A 5 ), 6 12 = A 3 , (55) 
the potential in Eq. ([3]) becomes V = Vq + xfFo, where Vb = nfxi + fi% x 2, and 

F = b n + 6 12 6 + 6 22 £ 2 2 + + ^4 , & = — , i = 2, 3, 4. (56) 
Using the Cauchy-Scwartz inequality 

2 



$1*2 



< $l$x$t$ 2 = i^i 2 !^! 2 , (57) 



we obtain the condition 

/(6,6,6) = 6 2 + £f-6<o. (58) 

We now introduce the Lagrange multiplier Ai in the quartic sector of the potential [l4[ related to the condition 
imposed by the Cauchy-Schwarz inequality Eq. (|58[) 

^(6,6,6^1) = F (6,6,6) + Ai/ (6, 6,6) = b n + (612 - A x ) 6 + + (633 + Ai)£ 2 + (6 44 + A x )^ (59) 

and apply the stability (positivity) condition in Eq. (|59[) after obtaining the derivatives 

OF dF dF OF „ 
<?6 96 »6 0A1 

to evaluate the minimum value for F (6, 6, £4, Ai) in the region of interest. 
The following equations are to be solved. 

26226 + h2 - A x = 0, 2 (fo 33 + Ax) 6 - 0, 2 (fe 44 + Ax) ^ = 0, $+£-6 = 0. (61) 

There are two solutions denoted by Ai, i = 1, 2, where 6 = £3 + £4, and 

Ai = 26226 +612, (633 + 26326 + 612)6 = 0, (644 + 26226 + 612)6 = 0, (62) 

where the stability condition becomes 

F (6, 6,6, Al)| min = (6ll - 6226" + (&33 + Al) £ 2 + (644 + Ai) 6 2 ) Imin > 0. (63) 



There is another solution, case B, in which Ai = 0, and £f + £f — £2 < 0. Where 

F (6,6, 6)l min = (611 + &126 + 6 22 6 2 + ^.3 + 6 44 6 2 ) Uin > 0. (64) 
In the first solution A\ we consider: 

6 = 0, 6^0, £ 2 = --L(& 12 + 6 33 ) = 6 2 - (65) 

^022 

The conditions and implications for a minimum in the region of interest are: 

612 + 6 33 < 0, 622 > 0, A 2 > 0, A3 + A 4 + A 5 < 0. (66) 

In the second solution A2 

6^0, (633 + 26226 + 612)6 = 0, 644 + 26226 + 612 = 0. (67) 
We have two possibilities, 6 = 0, and 6 0, and the existence of a minimum requires, if 6 = 

6=6", 6 = - {hl \t hi) , 622 > 0, 6 12 + 6 44 <0, A 2 >0, A 3 + A 4 -A 5 <0. (68) 
^022 

If 6 7^ 0, the implications are 

(612 + 644) (612 + 633) 
6 = t^t = t^t , (633 - 644) 6 =0, A 5 = 0, A 3 + A 4 < 0. (69) 

^t>22 ^022 

In case B, the equations to solve are 26226 + ^12 = 0, 26 33 6 = 0, 26446 = 0, and the conditions to have a 
minimum with its implications are 

61 9 

6 = -7^- > 0, 612 < 0, 6 22 > 0, A 3 < 0. (70) 

^022 



Now we apply the stability condition in Eq. (|63j) and obtain in cases A\ and A2 

4&22&11 > (612 + 6 33 ) 2 < (A 3 + A 4 + A 5 ) , (71) 

46ii6 22 > (612 + 6 44 ) 2 -2^X^ < (A 3 + A 4 - A 5 ) , (72) 

6 33 = 644, A 5 = =» -2VA1A2 < (A3 + A 4 ) . (73) 
Now, in case B and Eq. (|rJ4"|). the result is 

46 22 6n > (612) 2 =>> A 3 > -2y/\^. (74) 

Performing the second derivative 

F (6,6,6, Ai) = 6n + (612 - Ai) 6 + 6 22 6 2 + + Ai) 6 2 + (644 + AO g. (75) 

We obtain 

d fdF\ d 



, '(26226 + 612- Ai)= 26 22 >0. (76) 

36 \0i2J 06 

After imposing the stability condition, with one restriction we obtain these new boundary values for the quartic 
couplings: 



- 20^ < (A 3 + A 4 + A 5 ) , -2VA2A7<A 3 . (77) 
Considering now two restrictions: Following the same method as in the former case, considering now the conditions 

x\ + x\ < X1X2, xi+X2 — v 2 =0 (78) 
and two Lagrange multipliers, the function in consideration is: 
F 2c (xi,x 2 , x 3: xa, Ai, A 2 ) = 6nx 2 + 622^2 + ^33^3 + 644a; 2 + 612X1X2 + Ai (x\ + x\ - xix 2 ) + A 2 (x\ + x 2 - v 2 ) (79) 
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The new equations to solve are 



26iixi + 6122:2 + A 2 - Ai.T 2 = 0, 2622^2 + bi 2 xi + A 2 - A1X1 = 0, 

(6 33 +Ai)a;3 = 0, ( 6 44 + A x ) x A = 0. (80) 
Ai {x\ + x\ - x\x 2 ) = 0, A 2 (xi + x 2 - v 2 ) = 0. 

For F 2c (x 1) x 2 ,x 3) x 4 ,A 1 ,A 2 )\ min 

( 612 - Ai - 26 22 \ ^ n ( b 12 - A x - 2b u \ A 

x i = 77 — r 77 i-T? A 2 > 0, x 2 = — — A 2 > 0. (81) 

V4622611 - (612 - Ai) 2 / V4622611 - (612 - Ai) 2 ) 

and 

F 2c (x 1 ,X2,X3,Xi,A 1 ,A 2 )\ miD = bnxf + b 22 x?, + (b 12 - xix 2 \ m . n (82) 
The solutions satisfy the following requirements: 

F 2c (x 1 ,x 2 ,x 3 ,Xi, Ai, A 2 )| min _ (6n + 622 - 612 + Ai) 



with 



The requirements in Eqs. (1811) imply 



A 2 4fc 2 S„ -(S 12 - A,) 2 >0, (83) 



A; w,, -(»■;- A.)' , 3<0 . (84) 

2 (612 - Ai - 622 - on) 



(612 - Ai) 2 - 4611622 < 0, 612 - Ai - 26 22 < 0, 
(612 - Ai - 611 - 6 22 ) < 0, 612 - Ai - 26n < 0. 



(85) 



The aditional solutions that come from Eqs. ([80]) yield 

Ai = -633, Ai = -644, 633 = 644 • 

Then 



(A 3 + A4±|A 5 |) 2 <4AiA 2 , A 3 + A 4 ± |A 5 | < Ai +A 2 , A 3 + A 4 ± |A 5 | < 2A 2 , A 3 + A 4 ± |A 5 | < 2Ai 



and 



-2x/ATA2"< (A 3 + A 4 ±|A 5 |), A 3 + A 4 ±|A 5 | <2Ai, A 3 + A 4 ± |A 5 | < 2A 2 . 
Thus we obtain previous results plus 

Aj + A 2 > A 3 + A 4 ± |A 5 | , Ai+A 2 >A 3 + A 4 . (86) 

Bounds from extreme stability conditions 

Let us now determine the behavior of the quartic couplings in the case where the quartic Higgs potential V4 has its 
lowest possible value. In correspondence with Eq. ^ using the notation of Eq. (fS"4"f when x\ = vf, x 2 = v\ the V4 
can be simplified as follows 

V 4 = A1.T 2 + A2.T2 + XtXix 2 = (\/XiXi - \/A2X 2 ) + (\t + 2a/ AiA 2 ) xix 2 > 0. 



In the Extreme case, the condition to be satisfied is V4 — 0, then 
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In this case the Higgs masses become 




1 /4 

(4AiA 2 ) ' v, Ai 7- A 2 , v l ^ v 2 , 

(2A) 1/2 v, Xi = X 2 = A, wi = v 2 = v/y/2. 



In another interesting case, which is the Semi-extreme case, the V± = (\/AiXi — \/Ai£Ci) > 0, and 

^(tan/3) 4 , A T --2VAiA^, 

the M H o becomes 

M H o = V2 (Ax« 2 + (A 2 - Ai) vl) 1/2 • 

In both cases 

/I \l/2 

M fe o = 0, M ff± = ( - |A 4 + A 5 |j v, M A o = \X 5 \ 1/2 v. 



Numerical evaluation of the Higgs masses in terms of tan /3 



In general, according to Eqs. (f29)) -(f3"T j) the mass dependence on vi, « 2 , or v = y 'v\ + v%, can be reformulated in 
terms of v and u 2 . Therefore, with known w one can plot those masses in terms of u 2 , and determine their dependence 
on V2 or tan/3 (tan/3 = v-ij \/v 2 — uf), as in Fig. [TJ 
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FIG. 1. The V2 dependence of the Higgs masses for various sets of quartic couplings. 

As the favoured Standard Model Higgs mass window is still open, and as the 2HDM Higgses are nonstandard, 
almost none of the masses range has yet been ruled out, one can explore several hypothetical situations and analize 
the consequences of some published values for M H ± . 

We will now proceed to numerically evaluate the Higgs masses under different conditions for the quartic parameters 
at the energy scale E = M t , where M t is the mass of the quark top. First we will reproduce the M H ± value given 
in [l5j], considering several sets of A's. We choose A4 = —5.52, and A5 = —6.0 from many combinatios that give 
M H ± = 609.1, M^a = 621.7. As we shall see, in the extreme case, with one of the symmetries considered in [16[ in 
which Ai = A 2 = A, and tan/3 = 1, or tan/3 = 5 all the Higgs masses acquire constant values and M^o = 0. With 
the previous A4, A5, we obtain for H° values which are not ruled out experimentaly (for the SM Higgs) according 
to [l3, in the following way: (A = 0.485, A 3 = 10.55, A T = -0.97, M H o = 250), (A = 0.653, A 3 = 10.21, A T = -1.30, 
Mjja = 290) and (A = 1.94. A3 = 7.64, At = -3.88, M H o = 600). 

Though the masses do not depend on tan/3 at this energy scale, we will explore, in the next section, their behavior 
and dependence on it at higher energies scales. 



M A 

-M h , 

=179.47 
=226.06 
=234.45 



\ =12.5664, X =2.1608, X =1 1.6837 

1 '2 'a 

X =0.5104, X =-12.25 




v, [GeV] 
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Let us now classify the several cases, we will consider , in terms of the different stability condiions for the Ai s 
A.- Extreme case in which both equalities are satisfied 



A T = -2VA1A2, U A1/A2 = (tan/3) 4 



(87) 



Bl.- Semi-extreme case: 



B2.- Semi-extreme case: 



C- Lagrange inequality condition 



D.- Yukawa - Universality condition 



At 



-2v/A^, Ai/A 2 ^ (tan/3) 4 
A T ^-2VAiA^ A1/A2 = (tan/3) 4 
At > —2\/ Ai A2 



tan/3 = M t /M b , 



9t = 9b 



(89) 
(90) 
(91) 



The former cases will be combined with two additional conditions for the quartic couplings: 

Case 1 : Ai = A 2 , and Case 2: Ai ^ A 2 . 
In case 1A , 2A and IB, all the masses have a constant value in the inerval < u 2 < v due to their only dependence on 
v, but in spite of the an explicit independence of the masses on w 2 , tan /3 plays an important role on the energy scale 
dependence of the masses, as we shall see in the numerical solutions of the Renormalization Group Equations. which 
is fixed. In case 2B, M H o depends explicitly on w 2 and therefore on tan/3. 

To compare with the masses M H ± = 609, Ma — 621.7 given in Ref . [la] . we consider three kinds of compositions of 
the Ai i — 1, 5 parameters as in a), b) and c), which produce seven cases. The properties of these cases are to be 
analized in the following section. 

Ai A 2 A3 A4 A5 
0.485 0.485 10.55 -5.52 -6.0 



a. -When 
obtain 



which means that At 



-2y%Ai = -0.97 and (Ai/A 2 ) 



1/4 



1.0, we 



Case tan/3 M h o M H o w 2 
1A 1.0 250 179.47 
ID 41.2 250 253.7 



b.- When 



c- When 
tan /3 cases 



Ai A 2 A3 
6.0 1/6 9.51 



Ai A 2 A3 
12.0 3.0 11.5 



A 4 

-5.51 



A 4 
-5.5 



A 5 

-6.0 



A 5 

-6.0 



-2y%J^ = -2.0 and (Ai/A 2 ) 1/4 = 2.45, we obtain 



Case tan/3 M^o Mjjo u 2 
2A 2.45 358.9 234.9 



with A T = 0.0, = -12.0 and (Ai/A 2 



,1/4 



1.4, we explore different 



Case tan/3 M^o M^a w 2 

2B2 1.4 507.6 717.9 207.2 

2C 2.0 556.0 556.0 227.0 

2C 5.0 243.9 609.6 248.9 

2D 41.2 30.2 621.5 253.7 

d.- Considering now smaller values for {M H ±, Ma} = {253.8, 240.8}, with interesting properties for the energy range 
of validity of the 2HDM, with A T = -2^17^ = -0.97 and (Ai/A 2 ) 1/4 = 1.0, for ( ^ g Q A | 8 ^ ^ _ A Q 5 9 

Case tan /3 M^o Mjjo w 2 
1A 1.0 250 179.5 
ID 41.2 250 253.7 
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e.- Finally, for even lower masses {M h ±,Ma} = {170.3, 160.5}, which arise from ( q^25 q 125 0^65 

where At = — 2VA1A2 = —0.25 and (Ai/A 2 ) 1,/4 = 1.0, we obtain 

Case tan (3 M^a M H o V2 
1D,1B1 41.2 127 253.7 

VI. TRIVIALITY CONSTRAINS. 
Renormalization group equations 



A4 
-0.5 



A 5 

-0.4 



In this section we explore the asymptotic behavior of the parameters in the model, and their relations, through the 
Renormalization Group Equations (RGE) [8]. The RGE are a powerful tool to determine by the triviality principle, 
the energy bounds of the parameters and the validity of the model. In order to proceed in this way, to numerically 
evaluate the energy dependence of the A^ quartic couplings, it is necessary to consider the RGE of all the parameters, 
i.e., the gauge group couplings g±, 52, 33 of the symmetry groups U(l), SU (2), SU(3), the vacuum expectation values 
Vi, V2, and the Yukawa couplings of the top and the down quark sectors gt and ga respectively Refs. [HI]. 

The RGE determine the dependence of the coupling constants and other parameters of the Lagrangian on t, defined 
as t = In (E/m t ), where E is the renormalization point energy. The RGE for the gauge couplings gi, g2, <? 3 are: 



dgk 
dt 



1 



(4tt) : 



7 b k gl (i = 1,2,3), 



The RGE for the Yukawa couplings of the top and bottom quarks gt,gb are 

dg t 1 /9 



dt (4tt) 

dgb 1 
dt 



9 2 

(4tt) 2 \ 2 9b 



|9 b 2 -(^3i + ^52 2 + 853 2 ))'/ 



1 1 Q 

2 / 2 2 

TTfft - (t5i + 752 



% 3 2 : 



and for the vacuum expectation values v\ and V2 

d 1 



dt Vl 



di V2 



(Aw) 2 
1 

(4^)2 



[-3a 



((9/20) gj + (9/4) ,g 2 2 )] v x , 
-((9/20). g 2 + (9/4). 92 2 )] v 2 



(92) 

(93) 
(94) 

(95) 
(96) 



In the equations for the quartic couplings we include the quark Yukawa contributions of both sectors. 



d\x _ 1 
~dt ~16tt 2 
d\ 2 1 
dt ~16tt 2 



{24 (Ax) 2 
{24 (A 2 ) 2 



3Ai 
3A 2 



3 5 2 
3 5 2 



(.9') 



+ 2(A 3 ) Z + (A 4 )" + (A 5 )" + 2A 3 A 4 
+ 2(A 3 ) 2 + (A 4 ) 2 + (A 5 ) 2 + 2A 3 A 4 



{9') 



(9') 



^ =^2 { 4 ( A 3) 2 + 4 (3A 3 + A 4 ) (X 1 + A 2 ) - 3A 3 [3. 9 2 + [g'f - 2 ( 5 2 + g 2 b ) 



+ 2(A 4 ) 2 + 2(A 5 ) 2 + 



(gf 



^9t9 2 b 



(97) 



=7^2 { 4 ( A 4) 2 + 4A 4 (Ai + A 2 + 2A 3 ) - 3A 4 [s 9 2 + (g'f - 2 (g 2 + g 2 )] + 8 (A 5 ) 2 + 3.g 2 (g'f + 12g 2 g 2 } 



^ "ii A5 { 4 (Ai + A2 + 2Aa + 3A4) " 3 [ ig2 + {9 ' )2 ~ 2 ^ + 



9l) 



)]} 



The former equations are the coupled, non linear, ordinary differential equations whose solution provides the informa- 
tion about the renormalization point energy dependence of the masses of the five Higgs particles of the 2HDM. To nu- 
merically solve the RGE, the initial or final conditions for the parameters have to be previously chosen. In order to do so 
we use Ref. [2j. The range of values, we take, for the energy and the variable t are (Eq = M t ,E u ) = (173.2, 1.234 • 10 13 ) , 
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(to = 0, t u = 25) respectively, where M t stands for the mass of the quark top and E u corresponds to the electroweak 
unification energy where g\{E t ) — g 2 (E t ). The gauge couplings (gi, g2,gz) E=M — (0.4627,0.6466,1.2367,) are ob- 
tained using the following relations 

a e [M t ) = g 2 e /An = 1/ (127.9) , 9l (M t ) = y/5/Zg e / cos6 w 
g 2 (M t ) = g e / sm6 w , g 3 (M t ) = y/A-Ka s {M t ) 

where 9w is the Weinberg angle and sin 2 9w(M t ) — 0.235 and a s — 0.1217. The vev standard value that arises from 



v = 2M z /^gl+gl 



is v{M t ) = 253.81 GeV at M z = 91.19 GeV. 

In order to specify more rigorously the energy limits for the quartic couplings, we have numerically solved the RGE 
for the gauge group couplings gx, #2, 53, (Fig. [2]), the vacuum expectation values Vi, v 2 , and the top and the down 
quark Yukawa couplings gt and gd, under the following assumptions: 



The heaviest quark masses are related with the vevs v\ and v 2 and the Yukawa couplings and g^ 

V 2 n V 2 V\ 

M t = -j=gt, tan/3 = — , M b = -j=gb- 
V2 vi V2 



• The gauge bosons masses are related with the gauge couplings g' and g 

-v\j 'g 2 + {g'f \ 



M 1 AT MW 

M w = t; v 9i m z 



2 cos 9w 2 

where 9\y is the Weinberg angle and e the electron charge 

e = g sin 9w = g cos 9w- 
• Unification of the Yukawa couplings at E = M t or at E u , i.e., gi, = g t , and tan/3 = M t /Mb- 



(98) 

(99) 
(100) 




t [ln(E/m,)] 



FIG. 2. The energy dependence of the gauge couplings in the 2HDM. 



It is interesting to explore now, the energy bounds of the 2DHM, through the running of the quartic couplings which 
determine the mass values of the Higgses. In the case (c) considered in the previous section, when My± = 609, 
Ma = 621.7, the range of validity of the model is very short M t < E < 292 i.e., < t < 0.52 as can be seen 
in Figs. [SHHl There is an intermediate class of the models depicted at Figs. [71 [51 which have an intermadite range 
of validity < t < 11. So we will rather focus our atention on the cases where we can explore the universality 
of the Yukawa couplings and its unification, to study the mass-hierarchy problem. In this case, as can be seen in 
Figs. l9TfTT| the 2HDM is valid in the whole range of energies, this means Mt < E < E u where E u is the electroweak 
unification energy. In Fig. [3] we observe very slow dependence of the quaric couplings and the Higgs masses on the 
renormalization point energy. This model is characterized by rather small values of the quartic Yukawa couplings and 
the value of tan/3 such that it permits the unification of the Yukawa couplings of the up and down quarks gt = gb- 
In Figs. [TU1 HH we show the evolution of the Yukawa couplings, quartic couplings and the Higgs masses for the case 
when the Yukawa couplings are unified. In Fig.[TU]we assume that they are unified at low energy and in Fig. II II they 
are unified at high energy. The evolution of the quartic couplings and Higgs masses are similar in both cases. 
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t [InlE/m,)] t [ln(E/m t )] 

FIG. 3. The energy dependence of the quartic couplings and Higgs masses, case 2B2, with tan/3 = 1.41. 




t[ln(E/m ( )] tllnfE/m,)] 



FIG. 4. The energy dependence of the quartic couplings and the Higgs masses, case 2C with tan/3 = 2.0. 




t[ln(E/m ( )] tllnfE/m,)] 



FIG. 5. The energy dependence of the quartic couplings and the Higgs masses, case 2C with tan/3 = 5.0. 
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FIG. 6. The energy dependence of the quartic couplings and the Higgs masses, case 2D with tan/3 = 41.2. 




FIG. 7. The energy dependence of the quartic couplings and the Higgs masses, case 1A with tan = 1 (2D) case. 




FIG. 8. The energy dependence of the quartic couplings and the Higgs masses, case ID with tan/3 = 41.2 (2D) case. 




FIG. 9. The energy dependence of the quartic couplings and the Higgs masses, case 1B1 with tan/3 = 41.2. 




tllnlE/m,)! imam)] I [ln(EAnJ] 



FIG. 10. The energy dependence of the Yukawa couplings, quartic couplings and the Higgs masses in the tan/3 = 41.2 
case with Yukawa couplings g t = gb at low energy. 




FIG. 11. The energy dependence of the Yukawa couplings, quartic couplings and the Higgs masses in the tan /3 = 41.2 
case with equal Yukawa couplings g t = gb at high energy. 
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VII. RESULTS AND CONCLUSIONS 

With the aim to explore the Higgs mass content of the 2HDM extension of the standard model, among the different 
forms of the Lagrangian describing the same physical reality, we have chosen a specific one, in which the vacuum 
expectation values of both Higgs fields are real, and for simplicity also preserving the CP symmetry. We have deduced, 
in this model, the analytical expressions for the masses of the five predicted physical Higgs particles, and expressed 
the Higgs potential in terms of those masses, using Eqs. ([50)1 and (pTTj). We have also obtained, through the mass 
formulas, a set of constraints to be satisfied by the scalar parameters that determine the couplings and self-couplings 
of the Higgs fields introduced in the potential Eq. §5$ , and through the vacuum stability principle plus the Lagrange 
Multipliers method, and obtained additional conditions to be satisfied by those couplings. 

Ai >0, A 2 >0, 4AiA 2 > (A3 + A4 + A 5 ) 2 , (A 4 + A 5 )<0, A 5 < 0, A 4 <|A 5 |, (101) 

and 

Ai + A 2 > A3 + A4 + A5, A 3 + A4 + A 5 + 2VA1A2 > 0, A 3 + 2v/AiA 2 > 0, 4(AiA 2 )^A^. (102) 

We have also looked upon extreme and semiextreme conditions on the Higgs potential and gave a clasification of the 
different cases we analized under the RGE. 

As many authors base their calulations in symmetry conditions, such as Ai = A 2 and others in a phenomenologycal 
study of special events, it is important to analize the consequences of such assumptions and we tried at least partially 
address this problem. 

There is a batch of data to be analysed right now in search of some of the favored mass region, and all of it should be 
examined in the near future. The results of this paper may shed some light on physics of the Higgs sector depending 
on the properties of the Higgs particle. 

We have considered here, symmetries in the Ai parameters, universality of the Yukawa couplings at low energy Eq 
(M t scale) or high energy E u (weak-unification scale) , hierarchy of the quark masses and the energy range of validity 
of the model. The other symmetry considered here is the unification of the Yukawa couplings. It seems this symmetry 
makes the Higgs sector very stable as can be seen in Fig. HI 

In summary, the results in this paper may be a basis for further investigation in relation to the behavior and energy 
dependent characteristics of the Higgs particles. 

We finish by allegorically saying, that our paper still contains a "blank page" , which can only be filled after the 
discovery of Higgs bosons. 
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